Rigorous approach to the problem of ultraviolet divergencies in dilute Bose gases 
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In this Letter we consider a system of N pairwise finite-range interacting atoms and prove rigor- 
ously that in the zero-range interaction limit all the eigenstates and eigenenergies of the Hamiltonian 
converge to those corresponding to N atoms interacting via the Fermi-Huang regularized pseudo- 
potential. Next, we show that the latter eigensystem (if treated exactly) is invariant under a 
nontrivial transformation of the interaction potential. Finally, we realize that most of the approxi- 
mate schemes of many-body physics do not exhibit this invariance: we use this property to resolve 
all inconsistencies of the Hartree-Fock-Bogoliubov variational formalism known so far. 



The realization of Bose Einstein Condensate Q has 
brought an enormous interest in developing new theoret- 
ical approaches and refining the existing ones. The mean- 
field formalism with contact interactions has been shown 
to provide a powerful tool for analyzing the properties 
of trapped Bose gases Q . Unfortunately, the most gen- 
eral variational mean field approach: the Hartree-Fock- 
Bogoliubov approximation (HFB), is not yet quite satis- 
factory if used with contact interactions: it exhibits UV- 
divergencies, inconsistencies with the Hugenholtz-Pincs 
theorem ||, many-body T-matrix calculations and 
even with the very existence of atomic condensates them- 
selves ||. Several heuristic modifications of the theory 
were suggested HHHf^] , showing a good agreement with 
the experimental data J8|. In this Letter, we propose a 
novel, quite straightforward way to cure these inconsis- 
tencies. 

First, we prove rigorously that the regularized Fermi- 
Huang pseudo-potential ^ is not just an ansatz, but 
provides the exact zero-range limit of the many-body 
observables along with a cancellation of all the UV- 
divergencies. Second, we introduce a new family of 
pseudo-potentials parametrized by a free parameter A 
(so-called A-potentials) : no exact (after the zero-range 
approximation has been made) observable depends on it, 
while some approximate treatments differ for different A. 
The above conclusions are general and they do not rely 
on any particular approximation. Finally, as an appli- 
cation of this new potential, we find a particular value 
for A such that HFB equations are entirely free of all 
inconsistencies known so far. 

An explicit expression for the A-potential reads 
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V A (f) = g A S(f)[d r + A]{r-) with g A 



1 - aA 



(1) 



where f = r\ — r-i is the relative coordinate of two 
atoms, a is the s-wave scattering length, go = 2nh 2 a/jl 
is the usual effective coupling constant and jl = m/2 is 
the reduced mass. When A = 0, V A=0 coincides with 
the Fermi-Huang pseudo-potential. For a 1/r-divergent 
wave-function tjj(r) — d/r + rj(f) jlflL the action of the 
A-potential is 



V A iP(r) = g A S{r) r)(0) + Aa 



(2) 



For a low-energy two-body body scattering process, the 
eigenstates of the A-potential coincide with the asymp- 
totic form of the s-wave eigenstates of any other inter- 
action potential of a scattering length a. However, for 

h 2 

energies of the order - or higher, the actual finite size 

mar 

structure of the potentials comes into play, and the range 
of the applicability of the zero-range approximation re- 
duces to the particular case of a zero-energy resonance. 

Now, we consider a system of N atoms of mass m and 
coordinates {r} = {ri, ?2, • • • j Fn} interacting via V A . 
The Hamiltonian reads 
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where a — 1, .. . , a ma x = N(N — l)/2 labels ordered 
pairs of atoms (i a , j a ), and f a = Ti a — fj a is the relative 
position of the members of the a-th pair. As a direct 
consequence of Eq.(|2|), any eigenstate ^({r}) of (||) is a 
solution of an interaction-free Schrodinger equation sub- 
ject to the following contact conditions for all pair a: 
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where {7Z a } = ^R a , i ^ ia, ja} j is a set composed 

of the coordinate of the center of mass R a — (r ia +rj a )/2 
of the a-th pair and all other coordinates not belong- 
ing to this pair. Indeed, it is easy to show that the 5- 
singularities in the action of the Hamiltonian (||) on a 
many-body state ^({r}) cancel each other if and only if 
i$> satisfies the contact conditions (|J). Notice that these 
contact conditions do not depend on A, and thus no ex- 
act (after the zero-range approximation has been made) 
eigenstate does. 

Now, we are going to show how regularized pseudo- 
potentials arise in the limit of zero-range interactions. 
For this purpose, we consider N particles of mass m, 
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interacting via a potential which belongs to a one- 
parameter family of square- well potentials: v R (r) = 
—v R @(R—r). The depth of the potential v R > is 
chosen in such a way that the scattering length a is the 
same for all members of the family, and each of them 
supports the same number of s-wave bound states, ei- 
ther one or none depending on the sign of the scattering 
length . The Hamiltonian of the system reads 
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We wish to prove the following 

STATEMENT.- In the limit of infinitely small potential 
range, the Green's function of the finite-range-interaction 
Hamiltonian ^) converges to the Green's function of the 
pseudo-potential Hamiltonian (0): 



lim (E 
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Proof. As we saw above (see (Q) ) , all the A-potentials lead 
to the same eigenstates (and thus the same Green's func- 
tion): hence, without loss of generality, we can limit our 
proof to the case of the Fermi-Huang pseudo-potential 
V°. 

Let us define two operator- valued functions: 

Q E = (E + ie-a)- 1 . t E % = (1 - bQ? )~ l b . (7) 

The former is the retarded Green's function at energy E 
for a Hamiltonian a. The latter is the T-matrix of a per- 
turbation b in presence of the background Hamiltonian a. 
Two relations will be heavily used in what follows. First 
is the Lippman-Schwinger relation between the Green's 
function of the "full Hamiltonian" a + b and the one of 
the background: 



Q E - = Q E 

v a+b y " 
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The second is the Lupu-Sax formula [J12| relating the T- 
matrices of the same perturbation but in two different 
background Hamiltonians a\ and a.2'- 
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Introduce also a family of "reduced Hamiltonians" H. ? a 



and a family of reference Hamiltonians h E 
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where p a = (pi a — pj a )/2 is the relative momentum for 
the a-th pair. Each reference Hamiltonian is just a sum 
of the relative kinetic energy for the corresponding pair 
and the energy E at which the Green's functions (||) are 
compared. The Green's function of the a-th reference 



Hamiltonian is proportional to the zero-energy Green's 
function for the relative motion of two particles: 
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In turn the T-matrix of the interaction potential v R = 
v R (r a ) in presence of h E can be expressed through the 
zero-energy two-body T-matrix of it: 

({r} \t^ fiR \ {r'}) = flouVa.o 5{{n a n' a }). (12) 

The kernel D R is zero when r a > R or r' a > R and 
is normalized to unity as / d 3 fd 3 f ' D R (r, r') = 1. An 
explicit expression for it can be found in Ref. fll3|| . In 
the limit of zero-range interaction, the kernel obviously 
converges to a product of delta-functions, and hence the 
T-matrix converges to: 



goS(r a 
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Notice that by construction of the reference Hamiltonian 
h E neither the Green's function Qf„ nor the T-matrix 

T E E R depend on energy E. 

Using relations (||, ||), the full many-body Green's func- 
tion of the system can be rigorously expressed through 
the zero-energy two-body T-matrices ([l2]) of the inter- 
action potential v R (r). Removing from the Hamiltonian 
(||) one pair interaction after another, we obtain the fol- 
lowing chain relation: 
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\a} This relation is the cornerstone of the proof. 



Imagine now that the full Green's function G^, R 

acts on a state |\I>) whose wave function ^({r}) is 

regular everywhere. Computation of the result of 

this action ({r} involves expressions of a form 

%*2@*? a - 1} ~% ] %«S E Where S({r}) is a regular 
function, and in the limit R — > 0, leads to the expressions 
of the following type: 



r(*y - f 



(15) 
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(In Eq.(Ea), both T(r a ) and a depend also on {lZ a }.) T 
is given by 



r(ry= /d M {r'}({r}|^„ |{r'})c, <5(r a ')S({r'}), 

{a— 1} 



and a 



-go (/2/2^i 2 )S({r}) 



nition of the Green's function G^ R , we find that r(r Q ) 



r Q =0 



Now using the defi- 
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has an UV singularity of form a/r a (the same as the sec- 
ond term in the expression ( |l5| ) has) (see Jl^]). This leads 
to 



9oS(f a ) 



T(r a ) - ± 



(16) 



i.e. the expression (15) involves the Fermi-Huang 
pseudo-potential V° = V^{r a ) {c.f. (|) at A = 0). This 
justifies the following limit: 
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Inserting the above substitution at every level of the 
chain procedure (114) and collecting all the terms [[L5| one 



finally arrives at hni (?? R = Q^ p . p ., Q.E.D.. 

Notice that the relation (JTt]) clearly shows that the 
role of the regularizing operator in the pseudo-potential 
expression (|l|) is to subtract the free propagators Q? E 

already taken into account by the two-body T-matrix 
T® E As a result, UV divergencies disappear at each 

level of the chain recursion (|l4|). 

As an application of the A-potential, we consider now 
the HFB theory for N bosons interacting via V A with 
a > 0, in a box of a size L. As we will see, the A-freedom 
in choosing the effective Hamiltonian (||) offers the fol- 
lowing advantages: (a) Unlike for the conventional HFB 
formalism (A = 0), there exists a range of A such that 
the atomic condensate constitutes the minimum of the 
HFB functional in the low density regime [i6fl ; (b) For 
a particular value, A = A*, HFB equations are consis- 
tent with the results of the ladder approximation for the 
many-body T-matrix Q and the Hugenholtz-Pines theo- 
rem is satisfied; (c) In the vicinity of A* the ground state 
energy of the system is consistent with Bogoliubov's pre- 
dictions. 

The HFB approximation is twofold. First, it breaks 
the U(l) symmetry: the atomic field ip is split into a 
classical field $ and a quantum fluctuation <fi = V' — 
Second, the exact density operator is replaced by a Gaus- 
sian variational ansatz: D = exp(— K / 'k^T) / 'Z , where Z 
is the partition function and the quadratic variational 
Hamiltonian is 



K[h, A, 3] = - / / dWrJ [tf (fi)fe(ri,ra)0te) 

+ 0t(r 1 )A(r 1 ,r 2 )0t(r 2 )+h.c.l . (18) 



For what follows, we introduce the coordinates R = (fx + 
r 2 )/2 and r — r\ — r 2 . The second-quantized form of the 
full Hamiltonian (||) reads 



u A = ld 3 R ^t ( ___ A)v ,. 



where 



(19) 



F(R, R) = lim {d r + A}r F(R + r/2, R - r/2) 

r— >0 L 



(20) 



is a shortened notation for the action of the regular- 
izing operator (see ([[])). Using Wick's theorem, we 
obtain an approximate grand canonical potential J = 
E A - nN-TS, where E A = Tr[H A D] is the energy, TV is 
the number of particles, and S = — /«BTr[ln(Z))Z)] is the 
entropy Minimization of J with respect to the three vari- 
ational fields h, A and $ leads to the following implicit 
equations for these fields: 



h{n,T 2 ) = — (V 2 o)(r) + [KEi! n) 8{r) 
A(r*i,r 2 ) = hE A 2 8(f) 



(21) 
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A$ + [g A (2h + |$| 2 ) - f i] $ + .g A KA$* = , 



where KEii — 2ngA and hY^ A 2 = ,9a ( < I >2 + ^a) are the self- 
energies, n = Tr[0^ (i?)0(i?)Z)] is the non-condensed den- 

A 

sity, n = |<E>| 2 + n is the total density, and ka = k(-R, -R) 
results from the action of the regularizing operator (|2(i| ) 
on the anomalous density k(r\, r 2 ) = Tr[(/;(ri)</)(r 2 )Z3]. 

The diagonalization of the variational Hamiltonian K 
leads to the following quasi-particle spectrum: 
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Eqs.(gl|,g2j) clearly show that HFB is A-dependent. As 
it has been shown in Ref. 0, this approach is only able 
to provide a Born approximation for the diagonal self- 
energy TiEii, hence its explicit A-dependence. However, 
this is not the case for HT,±2', indeed, the total pairing 
field reproduces the contact conditions (f|) of a two-body 
wave function 

(V3(r 1 )v3(r 2 )) = ($ 2 + « )(l-^)+O(r) , (23) 

and as a result 7t£^ 2 = go(& 2 + «o) f° r a U A. 

Requiring that all the eigen-energies ( p2|) are real, we 
find that for zero temperature and densities below a value 

7T 

of n r ,it = : — o the existence of an atomic condensate 
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(4> =^ 0) implies the following constraint on A: 
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A*a < Aa < 1 with A*a 



«0 



(24) 



3> 2 + Kq 

At the lower limit A = A*, the A- regularized anomalous 
density disappears, and the theory become fully consis- 
tent with the results of the many-body T-matrix calcu- 
lations in the ladder diagrams approximation 
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(25) 



yielding a gapless spectrum JTt| , |18[ | . 

Consider now the zero-temperature low-density limit of 
our equations. Assuming Aa to be of the order of V na 3 
and neglecting all the terms of order na 3 or higher, the 
energy E A is independent of A and coincides with the 
well known Bogoliubov's result 
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15V? 

The A-potential based variational HFB model is there- 
fore consistent with the perturbative Bogoliubov's ap- 
proach. As the density increases the parameter A* in- 

7T 

creases as well, and at a critical density n cr i t = =- 

J 192a 3 

we find A*a = 1: the energy diverges and the mean field 
treatment breaks down. 

Note in conclusion that the A-invariance described 
in our Letter holds even if the constant A is replaced 
by an arbitrary field A(R). The generalization of our 
HFB theory to the case of the trapped gases is thus 
straightforward: One has simply to fix A as A(R) = 
A*(kq(R), <& 2 (J2)) according to ( p4] ) at every point R of 
the trap. 

As a extension of this work, we mention that using a 
procedure similar to the 3D case, it is possible to obtain 
the low-dimensional analogs of the A-potential: 



V A 



[P) ~ m log(gAi?) 
h 2 



Hp) i 



log{qAp)p-^- 



Vft(*) = --F 
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fj, Aam - 1 



8{z) 



1 d , 
1+ 2Ad~ z { 



where q = e c /2, C is the Euler's constant, R is the 2D 
effective hard disk radius and am is the ID scattering 
length 
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